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Abstract 

Each family of Gauss hypergeometric functions 

fn = 2-Fi(a + em, b + e 2 n; c + e^n; z), 

for fixed Ej = 0, ±1 (not all Ej equal to zero) satisfies a second order linear 
difference equation of the form 

Anfn-l + B„fn + Cnfn + l = 0. 

Because of symmetry relations and functional relations for the Gauss func- 
tions, many of the 26 cases (for different Sj values) can be transformed 
into each other. We give a set of basic equations from which all other 
equations can be obtained. For each basic equation, we study the exis- 
tence of minimal solutions and the character of f n (minimal or dominant) 
as n — ► ±00. A second independent solution is given in each basic case 
which is dominant when /„ is minimal and vice-versa. In this way, sat- 
isfactory pairs of linearly independent solutions for each of the 26 second 
order linear difference equations can be obtained. 



2000 Mathematics Subject Classification: 33C05, 39A11, 41A60, 65D20. 
Keywords & Phrases: Gauss hypergeometric functions, recursion relations, dif- 
ference equations, stability of recursion relations, numerical evaluation of special 
functions, asymptotic analysis. 
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1 Introduction 



The Gauss hypergeometric functions 

_ (a + ew, b + e 2 n \ 

where 

v 7 71=0 

satisfy three-term recurrence relations of the form 

Ai/n-l + Sti/n + Cn/n+1 = 0. (1.3) 

Examples are given in [1, p. 558]. 

In this paper we consider recursions with respect to n for the cases 

ej = 0,±1, j -1,2,3, e?+ £ 2 +e 2^ . (1.4) 

and we study the condition of these recurrences by obtaining the regions in the 
z-complex plane where a minimal solution exists. 

A solution /„ of the recurrence relation (1.3) is said to be minimal if there 
exists a linearly independent solution g n , of the same recurrence relation such 
that /n/jn^0asn^oo. In that case g n is called a dominant solution. When 
a recurrence admits a minimal solution (unique up to a constant factor), this 
solution should be included in any numerically satisfactory pair of solutions of 
the recurrence. Given a solution of the recurrence, it is crucial to know the 
character of the solution (minimal, dominant or none of them) in order to apply 
the recurrence relation in a numerically stable way. Indeed, if /„ is minimal as 
n — > +oo, forward recurrence (increasing n) is an ill conditioned process because 
small initial errors will generally dominate the recursive solution by introducing 
an initially small component of a dominant solution; backward recurrence is 
well conditioned in this case. The opposite situation takes place for dominant 
solutions. 

The problem of determining the z-values for which the /„ functions are mini- 
mal or dominant is considered in this paper. For the different recursion formulas 
we give a companion solution to /„ which, together with /„, form a numerically 
satisfactory pair of the corresponding three-term recurrence relation. 



2 Basic recursion relations and their solutions 

There are 26 recursion relations for these functions for all choices of Ej . However, 
we can use several functional relations in order to reduce our study to few basic 
recursions equations. First, we have the symmetry relation 
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In addition, the following relations can be used 



iM"' b ;z)=(l-z)- a 2F 1 ( a > * (2.2) 



c ' " ' \ c 



a > & . _"\ _ /1 „\c-a-6 17 /'c - a, c - 6 _ 



2 F 1 ^" c ;z)=(l-zy- a - b 2 F 1 ^ "> c " (2.4) 

See [1, p. 559]. These relations show that a great number of the 26 cases are 
equivalent or follow from each other. In fact we can find 5 basic forms that need 
to be studied, while the remaining cases follow from the relations in (2.1)-(2.4). 
Few cases for which Ej — 2 will be treated: only those which follow directly 
from the 26 cases for which \sj\ < 1, j = 1, 2, 3. 

When c = 0, — 1, —2, . . . the Gauss hypergeometric function in (1.2) is not 
defined. Also, when a or b assume non-positive integer values, the series in (1.2) 
terminates. In the following we will not distinguish about these special cases 
and the general cases, and we assume that all representations of the functions 
to be given are well defined, and we will not specify that special values of the 
parameters should be excluded in the results. 

In Table 1 we give an overview of all possible recursions, and indicate the 
five basic forms. Observe that we take recursion in positive direction equivalent 
with recursion in negative direction, however we will need to distinguish between 
both directions when studying the asymptotic behavior of the solutions; see 
Section 4. Apart from the notation in Table 1, we will also use the notation 
(sign(ei) sign(e 2 ) sign(e 2 )) when \ej\ < 1; for instance, the case k = 2 will be 
also represented as (+ + 0). 



2.1 Selection of a second solution 

Once we have reduced the number of basic recursions to be studied to 5, we will 
give for each basic form the coefficients A n , B n and C n of the recursion relation 
(1.3) and we will study the character of /„ as a solution of the corresponding 
recurrence relation. For this, we will need to find a second solution of the 
relation which forms a satisfactory pair of solutions together with /„ (i.e. a pair 
which includes the minimal solution when it exists). 

This second solution is chosen from several connection formulas between the 
hypergeometric functions. For example, the functions 

j-, fa, b \ i_ c fa — c+1, b — c+1 \ , . 

2-Fi I ' c ; z\ and z L C 2 F 1 I ^_ ^ ; z J (2.5) 

satisfy the same differential equation; sec [9, p. 112]. Both functions in (2.5) 
also satisfy the same difference equation, when a suitable normalization for the 
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Tabic 1: Only 5 basic forms remain 



k 


£\ 


£2 


£3 


type 


comments 


1 


1 


1 


1 


= 13 




use 


(2 


4) 




2 


1 


1 





basic 


form 










3 


1 


1 


-1 


basic 


form 










4 


1 





1 


= 13 




use 


(2 


3) 




5 


1 








basic 


form 










6 


1 





-1 


basic 


form 










7 


1 


-1 


1 


= 16 




use 


(2 


3) 




8 


1 


-1 





= 2 




use 


(2 


2) 




9 


1 


-1 


-1 


= 6 




use 


(2 


2) 




10 





1 


1 


= 13 




use 


(2 


2) 




11 





1 





= 5 




use 


(2 


1) 




12 





1 


-1 


= 6 




use 


(2 


1) 




13 








1 


basic 


form 










14 















void 






15 








-1 


= 13 




change sign 


in 13 


16 





-1 


1 


= 12 




change 


signs 


in 12 


17 





-1 





= 11 




use 


(2 


2) 




18 





-1 


-1 


= 15 




use 


(2 


2) 




19 


-1 


1 


1 


= 7 




use 


(2 


1) 




20 


-1 


1 





= 8 




use 


(2 


1) 




21 


-1 


1 


-1 


= 9 




use 


(2 


1) 




22 


-1 





1 


= 16 




use 


(2 


1) 




23 


-1 








= 17 




use 


(2 


1) 




24 


-1 





-1 


= 18 




use 


(2 


1) 




25 


-1 


-1 


1 


= 3 




change signs 


in 3 


26 


-1 


-1 





= 2 




use 


(2 


4) 




27 


-1 


-1 


-1 


EE 15 




use 


(2 


4) 
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second function is chosen. To obtain this normalization we use the relation 

(a, 6. \ T(c)T(c-a-b) ( a, b 
2Fl[ ~ > z] -T(c-a)T(c-b) ^{a + b-c+V 1 - 



(2.6) 



T(c)T(a + b-c) / c - ac -b 
+ T(a)T(b) {1 ~ Z) ^{c-a-b+V 1 -* 



(see see [1, Eq. 15.3.6] or [9, p. 113]). By replacing z by 1 — z and a + b — c + 1 
by c we can write this in the form 



Q z X ~ c iF x ( a ~ ° + 2- b c~ 0+1 '> z 



(2.7) 



where 

r(a + i-c)r(b + i-c) r(c - i) r(a + 1 - c) r(b + 1 - c) 

~ r(o + 6-c+l)r(l-c)' y_ T(a)T(b)T(l-c) 

(2.8) 

Another connection formula to be used is 

gg::ffg:::ii (-^-^.Cr--.;r^) 

r( c -i)r(6-c + i)r(i- a ) x _ _ a _ b A - a, 1 - 6 
r(h)r( c -a)r(i-c) — 2fl 2- c ' 

(2.9) 

where |phase(— z)\ < ir. This relation follows from [1, Eq. 15.3.7], after we 
change z — > 1/z and use the result (2.4). 

We can take one of the two terms in the right-hand sides of (2.6), (2.7) and 
(2.9) as a second solution of the equation that is satisfied by the function in the 
left-hand side. We need to verify if this second solution is linearly independent 
of the first solution. 

When selecting one of these candidates for the second solution, we can skip 
the gamma functions and other terms that are constant in the recursion. Also, 
gamma functions of the form T(a — n) will be replaced by using the relation (see 
[9, p. 74]) 

r(o-n) = - r. (2.10) 

v ; sin7rar(n+ 1 - a) v ; 
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2.2 Basic forms and their solutions 



(2.12) 



Now we summarize the main results that will be obtained for each basic recur- 
sion. We provide the coefficients of the recursion as well as a satisfactory pair of 
independent solutions, giving information on the character of the /„ functions. 
In later sections we will obtain the regions of existence of minimal solutions 
by means of Perron theorem and we will identify the minimal and a dominant 
solution by analyzing the asymptotic behaviour of well chosen solutions. 

2.2.1 Basic form k — 2 

The (+ + 0) recursion relation reads 

A 2 (a + n,b + n)y n _ 1 + B 2 (a + n,b + n)y n + C 2 (a + n,b + n)y n+1 = 0, (2.11) 
where 

A 2 (a,b) = (c- a)(c- b)(c- a- b- 1), 
B 2 (a, b) = (c - a - b){c{a + b - c) + c - 2ab+ 

+z[(a + b)(c - a - b) + 2ab + 1 - c]}, 
C 2 (a,b) = ab(c- a- b+ 1)(1 - z)2, 
with solutions given by 

fa + n, b + n 

Jn — 2-f 1 

_ T(a + n+ 1 - c)T(b + n + 1 - c) / a + n, b + n 
9n ~ r(a + fe + 2n-c+l) 21 \a + b + 2n-c+r l ~ Z 

(2.13) 

The second solution is taken from the first term in Eq. (2.7). 

As we will later show, when z < the recurrence has no minimal solutions, 
whereas in compact domains that do not contain points of (— oo,0], /„ is a 
dominant solution and g n is minimal. 

This case has applications for Jacobi polynomials. Wc have 

(-; a )(i±2)",« (-„;?-;,). P.U) 

A representation with +n at the a and b places follows from applying (2.4). In 
the interval of orthogonality — 1 < x < 1, we have z < 0, and if x e [—1,1] the 
recursion relation of the Jacobi polynomials can be used for computing these 
functions in forward direction. Only the usual rounding errors should be taken 
into account. 

Notice also how, in addition, information iis obtained regarding the recursion 
satisfied by 

, „ fa + n, b + n \ ,„,., 
U = 2FA c + 2n ;z\. (2.15) 
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By replacing z by 1 — z and c by a + b — c (sec also Eq. (2.10)) we see that 

_ T(c + 2n) / a + n, b + n \ 

9n - T(n+l-b + c)T(n+l-a + c) 2Fl { a + b-c ; 1 ~ Z ) (2 ' 16) 

is also a solution of the same recurrence, f n being minimal when z is in compact 
domains not containing points in the interval [1, +oo). 

2.2.2 Basic form k — 3 

The (+ H — ) recursion relation reads 



A 3 (a + n,b + n,c- n)y n -i + B 3 (a + n,b + n,c - n)y n + 
C 3 (a + n,b + n,c- n)y n+1 = 0, 



where 



(2.17) 



(2.18) 



A 3 (a, b, c) = -(a - c)(a - c - 1)(6 - 1 - c)(b - c)zll, 

B 3 {a, b, c) = c[aU + c 2 V + c 3 UV], 

ci = (1 - z)(b-c)(b- l)[a- l + z(b- c- 1)], 

c 2 = b(b+l- c)(l - z)(a + 6z - cz + 2z), 

c 3 = c — 26 — (a — 6)z, 

C 3 (a, 6, c) = a6c(c - 1) (1 - z) 3 V, 

U = z(a + b - c + l)(a + b - c + 2) + ab{\ - z), 

V = (1 — z)(l - a - & + a6) + z(a + b - c - l)(a + 6 - c - 2), 

with solutions given by 

'a + n, b + n 

fn = 2-Fl I ; Z 

c — n 



_ (-z) n r(a + 1 - c + 2n) T(b + 1 - c + 2n) (2-19) 
9n ~ r(a + n)r(& + n)r(2 - c + n)r(l - c + n) X 

a - c+\ + 2n, b-c+\ + 2n 
iF x | ; z 

2-c + n 

The second solution is obtained from the second term in Eq. (2.7) and applying 
Eq. (2.10). 

Later we will show that /„ is minimal as n — > +oo on compact domains 
inside the curve r = — 9 + 6\/3cos \Q with <6 < \~k (Figure 1) and that 
it is dominant in compact domains outside the curve. The opposite situation 
takes place for g n . 
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As n — > — oo, which corresponds to case k = 25 (i.e., (— — +)), the roles of 
/„ and g n are reversed, that is: f n is dominant (minimal) inside (outside) the 
curve represented in Figure 1. 

Similarly as happened for the case k = 2, we can also obtain the condition 
of the recursion for a recurrence outside our initial target of 26 cases. From Eq. 
(2.19) we see that the present case is related to the recursion for (ei£2£2) = 
(2 21). In addition, using Eqs. (2.1) and (2.2) we additionally see that case 
(2 2 1) is related to (2 - 1 1) and (-1 2 1). 

2.2.3 Basic form k — 5 

The (+ 0) recursion relation reads 

A 5 (a + ri)y n -i + B 5 (a + n)y n + C 5 (a + n)y n+1 = 0, (2.20) 

where 

A 5 (a) = (c- a), 

B 5 (a) = 2a-c-(a-b)z (2.21) 
C 5 (a) = a(z - 1), 

with solutions given by 

(2.22) 



where we have used the first term in the right-hand side of (2.9). 

As we will later show, in compact domains inside the circle \z — 1| = 1, /„ is 
dominant and g n is minimal, while the contrary takes place in compact domains 
outside the circle. 

2.2.4 Basic form k — 6 

The (+0—) recursion relation reads 

A 6 (a + n,c- n)y n -\ + B 6 (a + n,c- n)y n + C 6 (a + n,c- n)y n+1 = 0, (2.23) 
where 

Ae(a, c) = z(a — c)(a — c — l)(b — c)[a + z(b + 1 — c)], 
B 6 (a, c) = c[a(a - l)(c - 1) + a(a - I) (a + 36 - 4c + 2)z + 

{b - c)(b + 1 - c)(4a - c - l)z 2 -(a-b)(b- c)(b + 1 - c)z 3 ], 
C 6 (a,c) = -ac(c - l)[a -l + z(b- c)](l - z) 2 , 

(2.24) 
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fn — lF\ 



a + n, b 



(1 — a — n, 1— b \ 
2 -c - 



with solutions given by 



fn 



a + n, b 



c — n 



_ T(a+l-c + 2n)T(b+l-c + n) I a + n,b | 
gn ~r(a + b-c+l + 2n)T(l-c + n) 2 1 ^ + b c + x + 2n > 

(2.25) 

where g n is selected from the first term of Eq. (2.7). 

In this case /„ is never minimal. A satisfactory companion of g n is always 
h n = 9n — fn- The g n is minimal in compacts domains inside the region between 
the inner and outer curves in Figure 2. Contrary h n is minimal in compact 
domains inside the complementary region. As n — > — oo (case k = 22, that is 

(— +) or equivalently (0 h), that is, k — 16), the role which was played by 

h n for the (+ — ) case is now played by g n and the role of g n is played by /„. 

From Eq. (2.25), we observe that information can also be obtained for the 
recurrence (10 2) and then (by (2.1), (2.2), and (2.4)) also for the recurrences 
(212), (12 2), and (112). 

2.2.5 Basic form k = 13 

The (0 +) recursion relation reads 

Au(c + n)y n -i + B 13 (c + n)y n + C 13 (c + n)y n+1 = 0, (2.26) 

where 

A 13 (c)=c{c-l){z-l), 

Bia(c) = c[c - 1 - (2c - a - b - l)z], (2.27) 
Ci 3 (c) = (c- a)(c - b)z, 
with solutions given by 



fn — 2F1 



b ■ z 

c + n ' 



_ (-!)"(! -z) n T(c + n) fc-a + n, c-b + n 
9n ~ T(c-a-b+l + n) 21 ^ c-a-b+l + n ^ Z 



(2.28) 



where we have used the second term in the right-hand side of (2.6). 

As we will show next, /„ is minimal in compact domains inside the region 
Sftz < 1/2, where g n is dominant. Contrary, in compact domains inside the 
region 3?z > 1/2 /„ is dominant and g n is minimal. 

When n — > —00 (that is, for the case k = 15, corresponding to (0 0—)) /„ 
is a dominant solution both for < 1/2 and > 1/2. The satisfactory 
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companion solution for /„ (the minimal solution) is chosen differently for the 
cases < 1/2 and > 1/2 (sec §5.1). 



3 Domains for minimal and dominant solutions 

Perron's theorem (see [10, Appendix B]) gives in the case of finite limits the 
following results. Let a and j3 denote 

a:= lim ^, (3 := lim (3.1) 

Let ii and t 2 denote the zeros of the characteristic polynomial t 2 +at+[3 = 0. 
If \ti\ 7^ |*2 1, then the difference equation (1.3) has two linear independent 
solutions /„ and g n with the properties 

fn+l , 9n+l , lo r,\ 

— — ~ti, ~ t 2 . (3.2) 

Jn 9n 

If |* x | = \t 2 \, then 

limsup \y n \~ = \ti\ (3.3) 

n — ^oo 

for any non-trivial solution y n of (1.3). 

In the following subsections, we give for the five basic forms the domains 
in the z— plane where |ti| ^ \t2\- In these domains there is a true distinction 
between the two solutions of (1.3). If |ii| > \t 2 \ then the solution /„ that satisfies 
the relation in (3.2) is a maximal solution and g n is the minimal solution. On 
the curves where |ii| = \t 2 \ the two solutions are neither dominant nor minimal, 
and recursion in forward or backward direction is not unstable. 

For all five basic forms the ratios A n /C n and B n /C n of the difference equa- 
tion (1.3) tend to finite limits as |n| — ► oo. Interestingly, in all cases these limits 
are functions of z, and they are not depending on the parameters a, b or c. 



3.1 The domains for basic form k = 2 

The limits a and of (2.11) are 

The zeros of the characteristic polynomial are 

' 1 = (i-W' <2= (TtW- (3 ' 5) 

The equation = \t 2 \ holds when z < 0, otherwise \ti\ > \t 2 \. 
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0.1 




Figure 1: The curve |ti| = \ts\ for the basic form k — 3. 

3.2 The domains for basic form k = 3 

The limits a and (3 of (2.17) are 

8z 2 + 20z - 1 „ 16z 



/? = -- ^. (3.6) 



(1-z) 3 ' ' (I-*) 3 
The zeros of the characteristic polynomial are 

l-20z-8z 2 + (8z + l)l l-20z-8z 2 -(8z + l)I 
2(1^ ' h - 2(1^)3 ' ( 37 ) 

We write this in the form 

27 - 18w 2 - w 4 + 8w 3 _ 32(1 + w) 



h 16(1 -z) 3 (3 + w) 3 ' 

_ 27 - 18w 2 -w 4 - 8w 3 _ 32(1 - w) 

h ~ 16(1 - zf ~ (3-w) 3 ' 



(3.8) 



where w — y/8z~+T. To find the curve in the w— plane defined by \ti\ = jfol, we 
write w = re 10 . This gives the curve described by 

r = -9 + 6\/3cosf0, -±tt < 9 < |tt. (3.9) 

In Figure 1 we show this curve in the z— plane. In the domain interior to this 
curve we have |ii| > \t 2 \- 

3.3 The domains for basic form k = 5 

The limits a and (3 of (2.20) are 

z-2 1 

a=- , /?=- . 3.10 

1 — z 1 — z 
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Figure 2: The curve = \t 2 \ for the basic form k — 6. 



The zeros of the characteristic polynomial are 

ti = l, t 2 = -^—. (3.11) 

1 — 2 

The equation = \t 2 \ holds when \1 — z\ = 1, which defines a circle with centre 
z = 1 and radius 1. Inside the circle we have | i 2 1 > 

3.4 The domains for basic form k = 6 

The limits a and (3 of (2.23) are 

z 2 -6z + l 4z 

a = -^l-^< ^"(T^F- (3 ' 12) 



The zeros of the characteristic polynomial are 

Az 
(l~z) 



h = l, h = -- —. (3.13) 



To find the curve defined by = \t 2 \ 7 we write z = re l6 . This gives the curve 
described by 



r = 2 + cos6>±Vcos2fl + 4cos6l + 3, -tt < 6» < tt. (3.14) 

Both signs give a closed loop with common point — 1. In Figure 2 we show 
this curve in the z— plane. In the domain interior to the inner curve we have 
\ti\ > \h |; between the inner curve and the outer curve we have |ti| < \t 2 \, and 
outside the outer curve |ii| > \t 2 \. 
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3.5 The domains for basic form k = 13 



The limits a and {3 of (2.26) are 



2z-l n z- 1 ,„„„. 
a= , /3= . (3.15) 

z z 



The zeros of the characteristic polynomial are 

z-1 

h = 1, t 2 = . (3.16) 

z 

The equation = \t 2 \ holds when = \. When > | we have |ii| > |*2|. 

4 Asymptotics for minimal and dominant solu- 
tions 

In some cases we simply use the power series in (1.2), which provides an asymp- 
totic expansion for large c. In some other cases we can use connection formulas 
for transforming the Gauss function to the case for large c. We also use the 
integral representations 

* t rmh / v ' (1 - t)M[l - tz) ~" * ( "> 

where in the first integral 3?c > Wd > and in the second one Wd > 0. In the 
second integral the contour starts and terminates at t = and encircles the 
point t = 1 in the positive direction. The point t = 1/z should be outside the 
contour. The many-valued functions of the integrand assume their principal 
branches: the phase of (1 — tz) tends to zero when z — > 0, and the phases of t 
and (t — 1) arc zero at the point where the contour cuts the real positive axis 
(at the right of t = 1). To prove (4.2) integrate (if 3?c > > 0) along the 
interval (0, 1) with proper choices of the branches of the many-valued functions 
at the upper and lower sides of (0, 1), and obtain (4.1); see also [9, p. 111]. The 
integral in (4.1) is Euler's well-known standard representation. 

We use saddle-point methods (see [11]) to obtain asymptotic estimates of 
these integrals when one of more parameters are large. We omit details of 
the saddle-point analysis because we only need to obtain the main terms in 
the asymptotic estimates for identifying minimal and dominant solutions. We 
assume that z is fixed and properly inside the domains described in Section 3. 

For recent papers on uniform asymptotic expansions of hypergeometric func- 
tions, see [5], [6] and [7]. 
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4.1 Asymptotics for basic form k = 2 

For f n of (2.13) we use (4.2) and we obtain 

where 0(t) = lni — ln(t — 1) — ln(l — fcz). The saddle-points are found by putting 
(j)'(t) = 0, giving t± = ±l/y/z. If z £ (0, 1) we have t+ > 1 and the saddle-point 
contour starts at t = 0, turns around t = 1 through t + , and returns to t = 0. 
The dominant term in the asymptotic estimate is e^*+) = t™ (see (3.5)). This 
also holds for complex values of z. 

For <7„ of (2.13) we use (4.1) and we obtain 

_ T(a + n + 1 - c) T(b + n + 1 - c) 
ff " " ^/> + n+l-c)r(& + n) X (4 4) 

i 6 " 1 ^ - t)^ 6 -^! - t(l - z))- a e"^« dt, 



/ 



where ^>(i) = lni + ln(l — t) — ln(l — t(l — z)). The saddle-points are t+ = \ft%, 
t- = y/t{. When z £ (0, 1) we have t+ G (0, 1) and this point gives the dominant 
contribution. The dominant term in the asymptotic estimate is e^*+) = 
This also holds for complex values of z. 

When z is not a negative real number we infer that for g n the main term in 
the asymptotics is t% and that for /„ the main term is t\. We conclude that 
in compact domains that do not contain points of (— oo,0], g n of (2.13) is the 
minimal solution and /„ is a dominant solution of (2.11). 

4.2 Asymptotics for basic form k = 3 

For /„ of (2.19) we use (4.2) and we obtain 



fn 



r(c-n)T(l + b-c + 2n) r (1+) 



/ i 6 " 1 ^- If-^il - tz)- a dt, 
Jo 



2iriT{b + n) 

(4.5) 

where <f>(t) = lni-2 ln(t-l)-ln(l-iz). The saddle-points are t± = (l±to)/(4z), 
w = y/l + 8z. If z £ (0, 1) we have t+ > 1 and the saddle-point contour starts 
at t = 0, turns around t = 1 through i + , and returns to t = 0. The dominant 
term in the asymptotic estimate is e^*+) = (-l) n 2~ 2n q (see (3.7)). This also 
holds for complex values of z. 

For g n of (2.19) we again use (4.2) and we obtain 



(-z) n T(a + l- c + 2n) 

27U r(a + n)r(l - c + n) , . 

t b ~ c (i - l)- b (l - tz)-"^- 1 e n *W <fi, 

/o 



/' 

./o 
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where tp(t) = 21ni — \n(t — 1) — 21n(l — tz). The saddle-points are t± = 
(-1 ± w)/(2z). When z e (0,1) we have t + £ (1,2) and the saddle-point 
contour is similar as for /„. The dominant term in the asymptotic estimate is 
e^*+) = (-z)- ,l 2- 2n t™ (see (3.7)). This also holds for complex values of z. 

Taking into account the asymptotics of the gamma function in front of the 
integrals in (4.5) and (4.6), we infer that the dominant terms in the asymptotic 
estimate of /„ is t% ; and for g n it is t". We conclude that in compact domains 
interior to the curve of Figure 1, g n of (2.19) is a dominant solution and /„ is 
the minimal solution of (2.17). In compact domains exterior to this curve the 
roles of /„ and g n are interchanged. 

4.3 Asymptotics for basic form k = 5 

For f n of (2.22) we use (4.1) and we obtain 



r( c ) 



^rWF^I ^(i-tr^d-te)--*, (4-7) 

If |1 — z\ < 1 the main contributions to the integral come from t— values near 
the end-point 1. In that case the dominant term in the asymptotic estimate is 
*2 (see (3.11)). If |1 — z\ > 1 the main contributions to the integral come from 
t— values near the end-point 0. In that case /„ = 0(l/n) as n — > oo. 
For g n of (2.22) we use (2.3), which gives 

f\ -.M+b-c-» r ( a + n + 1 - c ) P fl + a-c + n, l-b z \ 
= r(a + n) 2Fl { 2-c 'T-^J- (4 - 8) 

We can use the same method as for /„, and conclude that in compact domains 
inside the disk \z — 1| < 1, g n of (2.22) is a minimal solution and /„ is a 
dominant solution of (2.20), respectively corresponding with t\ and ti of (3.11). 
In compact domains outside the disk the roles of /„ and g n are interchanged. 

4.4 Asymptotics for basic form k = 6 

For g n of (2.25) we use (2.1) and (4.1), and obtain 

* " r(?/„)ra- + c 2 + t) II - ^ - (1 - => ^l ^"*" , dt « 4 - 9 ' 

where 4>{t) = hit + ln(l — t). The point t = \ gives the main contribution, and 
the dominant term of the integral is Ar' a . 

Next we consider h n : — g n — f n . From (2.28), (2.7) and (4.1) it follows that 

_ z x - c + n r(c - n - l)T(a + 1 - c + 2n)r(2 -c + n) 

" — r( a + n)r(6)r(i- c + n)r(i-6) x (4 1Q) 



/ 

Jo 



t b - c (l - t)- b [l - t(l - z)]-^- 1 e n ^ dt, 
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where ip(t) = hit — 21n(l — tz). The end-point t = 1 gives the dominant 
contribution z n / (1 — z) 2n to the integral. 

Taking into account the contribution from the gamma functions in (4.9) and 
(4.10), it follows that g n corresponds with t\ of (3.13), for all z, and h n with t2. 

Because /„ = g n — h n , this function is never a minimal solution of (2.23). 
In compact domains interior to the inner curve and outside the outer curve (see 
Figure 2), g n is a dominant solution and h n the minimal solution; in compact 
domains between the inner and the outer curve the roles of g n and h n are 
interchanged. 



4.5 Asymptotics for basic form k = 13 

For f n of (2.28) we have the estimate /„ = l + C(l/n) for all z. For g n we apply 
(2.4), and obtain 

9n ~ z" Y{c-a-b+l + n) ^ [c a b + 1 + n ' 1 " Z ) ' (4U) 

We conclude that in compact domains of < |, /„ is the minimal solution 
and g n is a dominant solution. In compact domains of JJz > \ the roles of /„ 
and g n are interchanged. 



5 The cases k = 15, k = 16, A: = 25 

The cases fc = 15, k = 16, k = 25 in Table 1 are special because we refer to these 
as 'change signs in other cases', and we do not consider them as basic forms. 
The recursion relations for the cases k = 15, k — 16, k — 25 are the same as 
those for k = 13, k = 12, k = 3, respectively, when we recur backwards, that is 
to — oo. The zeros t\ and ti of the characteristic polynomial do not change when 
we change the recursion direction. In this section we give linearly independent 
satisfactory solutions for the cases k = 15, k = 16, k = 25. 



5.1 The case k = 15 

We verify if the solutions given for the recursion relation of the case k = 13 
can be used as satisfactory linear independent solutions for the same recursion 
relation used in backward direction. 

The asymptotics for of /„ of (2.28) for large positive n follows easily from 
the power series. For negative n we can use (2.7) and for the term with Q we 
use also (2.4). This gives 



(5.1) 



1— c+n/i „\c—n — a — b 77 fl ft, 1 b \ 

QZ {1 - Z) 2i H 2-c + n ' Z J 
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where P and Q are as in (2.8) with n replaced with —n. Again the power series 
can be used for obtaining the asymptotics, and it follows that for /_„ of (5.1) 
we have /_„ = 1 + 0(l/n) asrn +oo, when \z/{l — z)\ < 1; otherwise, the 
dominant term in the asymptotic behaviour is \z/(l — z)\ n . 

For the second solution g n of (2.28) we obtain, after using (2.4), 

(-l)"^~ c r( C -n) / z \ n ( 1-6, l-a \ 
g - n= T{c-a-b + l-n) \— z ) ^{c-a-b+l-n' 1 -*)- ^ 

We can verify as for the first F— function in (5.1) that the F— function in (5.2) 
has dominant term [(1 — z)/z] n when \z/(l — z)\ < 1; otherwise, the dominant 
term in the asymptotic behaviour is algebraic in n. Because of the extra factor 
[z/(l — z)] n in front of the F— function in (5.2), it follows that /_„ and <?_„ 
are minimal or dominant in the same z— domain. Hence the {/n,<?n} of (2.28) 
does not constitute a satisfactory pair of linearly independent solutions for the 
recursion (2.27) in backward direction. As a consequence, we cannot use case 
k = 13 in backward direction for the present case k = 15. 

The second solution g n of (2.28) has been selected by considering the second 
term in (2.6). When instead we consider the first term, using (2.10), with c 
replaced with c — n, we can take as another solution of (2.27) the function h n , 
where 

= r(n + 1 - c + a)r(n + 1 - c + b) ( a, b 

~ n ~ T(n + 1 - c)r(n +l-c + a + bf 1 \n+ 1 - c + a + b ] 

(5.3) 

The F— function is 1 + 0(l/n) for large positive n in compact z— domains. 
Hence, h n can be used as a proper second solution together with /„ for the 
recursion relation of the case k = 13 in backward direction when \z/{l — z)\ > 1. 
We still need another function for the complementary domain. We consider for 
this the connection formula (2.9) use the first term in the right-hand side, and 
(2.10), and we take as a new solution of (2.27) the function j n , where 

T(a - c + n + l)T(b -c + n+1) 

j-n = (-1)™ x 

T(2 -c + n)Y{\ -c + n) 



(5.4) 



1-zJ V 2-c+n 



The F— function is 1 + 0(l/n) for large positive n in compact z— domains. 
Hence, j n can be used as a proper second solution together with f n for the 
recursion relation of the case k = 13 in backward direction when \z/(l — z)\ < 1. 

Summarizing, for the present case k = 15 we can use the recursion (2.27) of 
basic form k = 13 in backward direction with /„ of (2.28) and h n of (5.3) for 
z— values in compact domains of the half-plane > ^ where /„ is a dominant 
and h n is the minimal solution. In compact domains of < \ we can use /„ 
together with j n of (5.4), with /„ as a dominant and j n as the minimal solution. 
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5.2 The case k = 16 



We verify if the solutions given for the recursion relation of the case k = 6 can 
be used as linear independent solutions for the same recursion relation used in 
backward direction. 

The function /„ of the case k — 6 given in (2.25) is 1 + 0(l/n) asm — oo. 
This easily follows when we use the integral representation in (4.1). Using this 
integral g n of (2.25) it follows that for g n the dominant term is ij n , where i 2 is 
given in (3.13). 

It follows that for the present case k — 16 we can use the recursion of basic 
form k = 6, that is, (2.23), in backward direction with /„ and g n of (2.25). Inside 
compact domains interior to the inner curve (around the origin) of Figure 2 and 
exterior to the outer curve /„ is a dominant and g n is the minimal solution. In 
compact domains between the two curves the roles of /„ and g n are interchanged. 

5.3 The case k = 25 

We verify if the solutions given for the recursion relation of the case k = 3 can 
be used as linear independent solutions for the same recursion relation used in 
backward direction. 

When we apply (2.4) to both /„ and g n of (2.19), we see that the asymptotics 
for n — ► — oo follows from that of n — ► +oo, when the roles of /„ and g n are 
interchanged. 

We conclude that in compact domains interior to the curve of Figure 1, f n 
is a dominant solution and g n is the minimal solution of of (2.17) when used in 
backward direction. In compact domains exterior to this curve the roles of /„ 
and g n are interchanged. 



6 Numerical examples 

The power series (1.2) is very useful for numerical computations for z properly 
inside the unit disk. Transformations and connection formulas as in (2.2), (2.3), 
(2.6) and (2.9) can be used to cover large parts of the complex z— plane. Other 
connection formulas are available with other z— arguments and in fact for the 
computation of the Gauss function we can use power series with powers of 



1 



1 

z 



1 



1 



\-z z-1 



(6.1) 



For numerical computations we need convergence conditions like 



\z\ < p, |1 - z\ < p, 



< P, 



1 



< P, 



1 



1 



< P, 



1 



< P, 



(6.2) 



with < p < 1. 

Not all points in the z— plane satisfy one of these inequalities for a given 
number p. In Figure 3 we take p=\- In the dark area at least one of the above 
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2 


p = 0.75 







Figure 3: In the light domains around the points e ±7Tt ^ 3 none of the inequalities of 
(6.2) is satisfied. 



inequalities is satisfied. In the light areas, "around" the points e ±nl ^ 3 , none of 
these inequalities is satisfied. By choosing p closer to unity these light domains 
become smaller. 

For certain combinations of the parameters a, b and c, the connection for- 
mulas become numerically unstable. For example, if c = a + b, the relation in 
(2.6) is well-defined, although two gamma functions are infinite. By using a 
limiting procedure the value of 2^1(0, b;a + b; z) can be found. For c close to 
a + b numerical instabilities occur when using (2.6). See [2] for many examples 
and details. 

Other instabilities in the evaluation of the power series (1.2) may arise for 
large values of a and b. 

In [10, p. 71] an example is given how to compute a Gauss function with 
argument z — e"^ 3 , the point that is excluded from the convergence domains 
shown in Figure 3. Wimp considers the computation of 



I 1 

4 ' 
3 



,7Ti/3 



27re w6r(i) 
9[r(f)] 2 ' 



by using a Miller algorithm for the hypergeometric functions 

, (n + a, n + b \ 

fn = 2 F 1 l 2n + c ;z\, n = 0,l,2,... 



(6.3) 



(6.4) 



This recursion type was not initially included in the group of 26 discussed in the 
present paper. However as, discussed in §2.2.1, this recurrence can be related 
to the case k = 2 to conclude that /„ is minimal. Therefore, Miller's algorithm 
can be applied when a sum rule is provided, as done in [10, p. 71]. 
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We can also use basic form k — 13, with pure c— recursion. From §3.5 and 
§4.5 it follows that for the point z = e 7 ™/ 3 = \ + the solutions of the 

recursion relation (2.26) are neither dominant nor minimal. We use backward 
recursion for 2-Fi(d, b;c + n; z) with two starting values for n — 29 and n = 30. 
With these large values of c+n the power series converges fast. In 15D arithmetic 
we have computed the value of (6.3) with a relative error 2 x 10~ 14 . The exact 
value is 

0.883319375142724... + 0.509984679019064.. A . 
With recursion we obtain 

0.883319375142719+ 0.509984679019039* . 

In [10, p. 72] another example of the Miller algorithm is discussed for the 
basic form k — 13 (pure c— recursion). 

Several algorithms based on recursion relations for special cases of the Gauss 
functions have been published, in particular for computing Legendre functions. 
For recent papers, see [3], [4] and [8]. 
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